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ABSTRACT 


A  model  for  the  statistics  of  sea  clutter  has  been  developed  from  scattering  theory  and 
the  composite  surface-scattering  model.  The  model  postulates  that  sea  clutter  is  exponentially 
(Rayleigh  envelope)  distributed  for  gla  md  should  tend  toward  the  lognormal  distribu¬ 

tion  (in  particular  for  horizontal  polarization)  with  increasing  roughness.  The  lognormaiity 
of  sea  clutter  arises  from  the  tilting  of  the  slightly  rough  “patches”  by  the  large-scale  rough¬ 
ness  (undulating  surface). 

An  empirical  identification  of  the  statistics  of  sea  clutter  taken  with  the  Four-Frequency 
Radar  system  shows  that  in  general  the  distribution  of  sea  clutter  is  intermediate  between 
the  exponential  (Rayleigh  envelope)  and  the  lognormal  distribution.  However,  for  calm  seas 
and  small  sample  sizes  (less  than  about  200  independent  samples)  the  distribution  of  sea 
clutter  may  be  approximated  by  either  the  exponential  or  the  lognormal  distribution. 

The  first  five  central  moments  of  sea  clutter  (in  decibels)  have  been  calculated  for 
moderate  and  rough  sea  conditions. 


PROBLEM  STATUS 

This  is  a  final  report  on  one  phase  of  the  problem;  work  on  other  phases  is  continuing. 


AUTHORIZATION 


NRL  Problem  R02-37 
Project  A310310B/652A/A  R02101-002 
(formerly  A310310B/652A/1  R008-01-020) 


Manuscript  submitted  September  10,  1971. 


ON  THE  STATISTICS  OF  SEA  CLUTTER 


INTRODUCTION 

Radar  detection,  which  is  a  binary  detection  problem  (i.e.,  the  selection  of  one  outcome 
from  two  events),  has  developed  quite  rapidly  since  the  early  1940’s  and  its  progress  can  be 
followed  in  Middleton  (1).  However,  modern  radar  detection  theory  is  not  considered  to 
have  started  until  the  Marcum  (2)  and  Swerling  (3)  contributions,  and  presently  it  has 
reached  a  high  level  of  maturity.  Nevertheless,  most  formal  developments  still  only  treat 
explicitly  the  detection  of  signals  (deterministic  or  random)  strictly  in  Gaussian  noise  (4). 

Unfortunately,  many  of  the  geophysical  noises  encountered  in  practice  are  non-Gaussian 
(e.g.,  atmospheric  disturbances,  terrain,  and  sea  clutter),  and  to  implement  realistic  detection 
schemes  the  true  statistics  of  these  geophysical  noises  must  be  identified  a  priori. 

Recently  it  has  been  suggested  that  the  statistics  of  sea  clutter,  in  particular  for  high- 
resolution  radars  and  toward  grazing  incidence,  cannot  be  expressed  as  a  Rayleigh  distribu¬ 
tion  (exponential  in  power)  but  can  be  approximated  by  other  distributions,  among  them  the 
lognormal  distribution  (5). 

The  main  purpose  of  this  investigation  is  to  identify  the  statistical  properties  of  sea 
clutter,  both  by  electromagnetic  scattering  theory  via  the  composite  surface-scattering 
model  (6,  7)  and  empirically  by  statistical  analysis  of  sea  clutter  taken  with  the  Four- 
Frequency  Radar  (4FR)  system  (8),  which  transmits  at  428  MHz  (P  band  or  UHF),  1228 
MHz  (L  band),  4455  MHz  (C  band),  and  8910  MHz  (X  band).  The  composite  surface¬ 
scattering  model  already  has  been  quite  successful  in  predicting  the  mean  value  of  sea 
clutter  (9-11)  and  explaining  the  width  of  the  doppler  spectrum  of  radar  sea  echo  (12).  The 
potential  of  this  scattering  model  to  derive  other  statistical  information  of  sea  clutter  is  to 
be  explored  here. 

Anticipating  some  of  the  results  to  be  obtained  later,  it  is  possible  to  say  that  the  tilting 
of  the  slightly  rough  “patches”  by  the  large-scale  roughness  (i.e.,  the  undulating  surface)  is 
the  mechanism  that  generates  the  non-Rayleigh  statistics  observed  in  sea  clutter,  in  particular 
for  horizontal  polarization.  Also  in  the  empirical  identification  of  sea  clutter  taken  with  the 
4FR  system  by  the  Kolmogorov-Smimov  test  on  the  cumulative  distribution  and  by  the 
computation  of  the  first  five  central  moments,  the  distribution  of  sea  clutter  (for  a  radar 
pulsewidth  of  0.5  psecond)  is  found  to  be  intermediate  between  the  exponential  (Rayleigh 
envelope)  and  the  lognormal  distribution,  and  these  distributions  may  serve  as  the  limiting 
distributions  of  sea  clutter. 
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THEORETICAL  CONSIDERATIONS 
The  Radar  Detection  Problem 

In  binary  detection  problems  the  Bayes  criterion  is  used  (1,  4)  when  the  a  priori  prob¬ 
abilities  and  cost  functions  of  making  a  decision  are  known  (i.e.,  the  risk  in  making  the  wrong 
decisions  is  minimized).  The  implementation  of  the  Bayes  criterion  in  detection  leads  to  a 
likelihood  ratio  test 


A(V)=|(^“)>X*i,/fli8true* 


da) 


and 


A(V0=f(^)<X’if//o“  true, 


(lb) 


where  p(  V/Hx )  and  p(  V/H0 )  are  the  conditional  probability  densities  given  that  event  Hl  is 
true  (signal  plus  noise  are  present)  and  given  that  event  H0  is  true  (noise  alone  is  present), 
respectively.  The  threshold  level  X  is  a  function  of  the  cost  functions  and  the  a  priori  prob¬ 
abilities  of  the  sources. 

Since,  in  practice,  the  cost  functions  and  the  a  priori  probabilities  are  not  known,  other 
criteria  must  be  considered.  In  radar  detection,  the  Neyman-Pearson  criterion  is  most  fre¬ 
quently  used,  the  probability  of  false  alarm  PF  is  specified,  and  the  probability  of  detection 
PD  is  maximized. 

The  threshold  level  under  these  constraints  is  obtained  from 


‘TO 

PFm  {  p(V/H0)dV, 


(2) 


and  the  required  signal-io-nois*  ratio  required  for  detection  then  is  obt'ined  from 


PD-  I  p(V/Hx)dV.  (3) 

\ 

The  implementation  of  the  Neyman-Pearson  criterion  also  leads  to  a  likelihood  ratio  test 
similar  to  Eq.  (1)  with  the  threshold  level  now  being  determined  from  Eq.  (2). 

Thus,  for  radar  systems  operating  in  the  sea  environment  the  conditional  probability 
density  of  clutter  must  be  known  in  order  to  have  realistic  climates  of  the  probability  of 
false  alarm  as  a  function  of  radar  and  sea  parameters.  Othe  c;  iteria  for  detection  also 
exist  (1,  4),  but  they  will  not  be  diacussed  here. 
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Integration  of  several  radar  pulses  should  improve  the  performance  of  deter  tion,  the 
amount  of  improvement  depending  on  the  correlation  properties  of  clutter  alone  m  relation 
to  the  correlation  properties  of  signal  plus  clutter. 

Electromagnetic  Scattering  from  Rough  Surfaces 

Consider  the  scattering  properties  of  a  statistical  rough  surface,  with  surface  currents 
K(  *  n  X  H  and  Km  *  -ti  X  E  due  to  a  plane  wave  incident  on  the  surface  (Fig.  1).  The 
scattered  fields  may  be  obtained  by  means  of  the  following  integral  equation  formulation 
(13, 14),  where  the  time  dependence  is  taken  to  be  eiu}t  and  the  unit  normal  vector  n  is 
pointing  toward  free  space: 


E,(r)-  ids'  Km(r')  XV rG(r/) v.y  j da'  K.(r')  X  V  'G(r,r') 


(4) 


H,(r)»-  f  dt  K^fr')  X  V'C(r,r’)  +  r~—V  X  f  da’ K„,(r')  X  V’G(r,r’), 
J.  0  1 


(5) 


where  to  is  the  electromagnetic  frequency  in  radians,  is  the  magnetic  permeability,  k  is 
the  propagation  constant  of  a  plane  wave  in  free  space,  G(r,r’)  is  the  Green’s  function  in  fr» 
space,  and  r’  is  the  radius  vector  Tor  a  point  on  the  rough  surface. 
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Pig.  1— pomtry 


4 


G.R.  VALENZUELA  AND  M.B  LAING 


The  integral  equations  (4)  and  (5)  can  be  simplified  for  the  “far  field”  (i.e.,  all  scattered 
rays  arriving  in  a  given  direction  must  be  propagating  parallel  to  each  other).  In  this  case, 
the  following  approximations  apply 


V'G(r.r')  *  V' 


(e-ik\r~  r'| 

\  4*r|r  -71 


)  30  ikG(rjr')n2 


(6) 


G(R,  r') 


p-ikR 

m~4wR 


eik2"r' 


O' 


VX  K4(r')  X  V  'G(r.r')  -  -  V  X  V  'G(r.r')  X  K.(r')  *  -*2 n2  X  (n2  X  K, )G(r,r'),  (8) 

where  k2  is  the  propagation  vector  of  the  scattered  wave,  R  is  the  distance  from  the  obser¬ 
vation  point  to  the  origin  of  the  coordinates  system,  and  n2  -  k2/k. 

Using  the  approximations  given  in  Eqs.  (4)  and  (6),  the  following  equations  result: 

Ef(i*)  *  K  j Jdt'  n2  X  lKm(r')  +  Z0n2  X  Ke(r')]e«‘2-'!  J  (9) 

H,(R)*-K{/ds'n2  X  [Ke(r')-y0n2  X  Km(r')]e*2*' j  ,  {10) 

where 

______ 

K  =  — -  ,  z0  "v/Mo/fo  *  120ft,  and  Y0  *  1/Z0. 

Accordingly,  the  statistics  of  the  scattered  fields  depend  only  on  the  statistics  of  the  surface 
currents  and  the  statistic  properties  of  the  rough  surface. 


Letting 


the  magnitude  square 


L  *  fds  Km(rVh2V 
M  -/da'  Ke(rV'2'f', 

of  the  fields  (which  is  proportional  to  power)  is  given  by 

2  f 


IL*(p  X  n2)|a  +  Z$|M-p|2  +  2Z0Se(M-p)(L*-(p  X  n2)] 


(11) 

(12) 


(13) 


and 


IE,-(n2  Xp)|2  *  2  JlL’pl2  +Z«|M-(n2  Xp)|*  ♦  2Z0 *e(L*-p)lM-(na  X  p)J  J  ,(14) 

where  L*  is  the  complex  conjugate  vector  of  L  and  p  is  a  polarisation  unit  vector  such  that 
p-n2  -  0. 


NRL  REPORT  7349 


5 


For  a  rough  surface  of  impedance  Zv  the  magnetic  surface  current  Km  (r')  can  be  ob¬ 
tained  from  the  electric  surface  current  Ke(r')  by 


Km(r')  =  ~Ztn(r’)  X  Ke(r'). 


(15) 


The  surface  currents  Ke  and  Km  are  determined  from  the  boundary  value  problem  of 
the  air/rough  surface  interface  from  the  incident  electromagnetic  wave  on  the  surface.  When 
the  scattering  surface  is  statistically  rough,  the  boundary  value  problem  cannot  be  solved  in 
closed  form,  and  accordingly  appropriate  approximations  must  be  used  for  the  surface 
currents. 

Fortunately,  for  rough  surfaces  like  the  sea,  which  has  a  two-dimensional  spectrum 
JF(K;  of  the  form  K~4  ( K  is  the  wavenumber  of  the  surface  roughness  or  waves),  a  slightly 
rough  local  boundary  condition  may  be  used.  11114  fundamentally  separates  the  sea  surface 
into  large-scale  roughnesses  (undulating  surface)  and  small-scale  roughnesses  (Bragg  resonant 
scatterers).  According  to  the  local -slightly  rough  boundary  condition,  the  surface  fields  on 
the  undulating  surface  (i.e.,  f(r’))  may  be  expressed  as  a  function  of  the  small-scale  rough¬ 
ness  (i.e.,  £(r')), 


H(r')-h[n(i'),k1^f)HrV*iV. 


(16) 


where  now  n(r')  becomes  the  unit  normal  vector  to  f(r'),  the  undulating  surface.  The  Soviets 
have  used  a  similar  approach  (15). 


Substituting  Eqs.  (15)  and  (16)  in  Eq.  (9)  and  specializing  for  backscatter,  k2  =  ~k1 
and  n2  =  -Hj ,  the  following  approximation  is  derived: 


Ef  «*X  na  X  /[nx  +  Z,n(r')]  X  (n(r')  X  Nr')tt(rV2,ki 


(17) 


The  integration  in  Eq.  (17)  can  be  performed  in  the  following  sequence:  first,  integrate 
over  each  slightly  rough  surface  “patch”  and  then  add  the  contributions  from  the  “patches” 


Ef  «  K  ^  Jnx  X  {«!  +  Z'ttfp))  X  (n(p)  X  h(p)]e-2*i*0>)f(p) 

P«teh«(p) 

[  £(r')e-2,(k  lx'x  +*iy'y  )  d(patch) 

(p) 


(18) 


The  backscattered  power  Pf ,  which  is  proportional  to  IE®  ia,  is  given  by 
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pf  {nj  X  [nx  +ztn{p)]  X  [n(p)  X  h(p)]  J  {ni  X  [rij  +  Zt n(£)]  X  [n(C)  X  h(£)]  • 

(P)  fC) 

g-2iifeix'(p)x  -*ixn(«)*,'+feiy'(p)y'-fciy’,(»)y'']  d(p)d(t),  (19) 

where  x'y  are  in  patch  (p)  and  x",y"  are  in  patch  (£). 

Neglecting  the  cross  ^rms  (the  contribution  from  two  different  patches),  basically  as¬ 
sumes  that  the  fields  from  different  patches  are  independent  and  the  backscatter  power  will 
be  given  by 

Pf  «  jlnl  X  Inl  +  Z.n(P)]  X  Jn(P)  X  h(P)l  I2' 

(P)  1 

/ /  t(r')£*(r'  +  Ar>"*l*i*'A*’+*iy,Ay,J  d(Ax')d(Ay')} • 
(p) 

(20) 

Thus,  for  an  incident  plane  wave  when  the  slightly  rough-local  boundary  condition  is  used, 
the  backscattered  power  from  a  rough  surface  is  given  by  Eq.  (19).  If  the  contributions  from 
the  various  “patches”  are  mutually  independent,  the  backscattered  power  may  be  approxi¬ 
mated  by  Eq.  (20),  which  can  be  given  the  same  physical  interpretation  as  the  composite 
surface-scattering  model  formulation  as  proposed  by  Wright  (6).  This  last  result  will  be  the 
basis  for  the  statistical  model  of  sea  clutter  to  h  developed  in  the  following  section. 

Here  it  is  appropriate  to  indicate  that  this  analysis  applies  for  an  incident  plane  wave.  In 
practical  applications  the  incident  electromagnetic  radiation  originates  from  an  aperture  of 
finite  dimensions,  and  the  incident  electromagnetic  radiation  is  a  superposition  of  many 
plane  waves.  Accordingly,  the  backscattered  power  should  include  an  additional  summation 
over  the  spectrum  of  the  incident  plane  waves. 

STATISTICAL  MODEL  OF  SEA  CLUTTER 

As  shown  previously,  the  scattering  results  for  certain  types  of  rough  surfaces,  in 
particular  for  the  sea,  car.  be  simplified  considerably  by  the  assumptions  of  the  com¬ 
posite  surface-scattering  model.  It  is  recalled  that  the  assumption  of  incoherent  addi¬ 
tion  has  been  used  (power  adds)  in  contrast  with  coherent  addition  in  which  the  fields  from 
the  various  slightly  rough  “patches”  add. 
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The  backscattered  power  from  each  “patch”  is  the  product  of  two  factors:  one  factor 
is  purely  a  function  of  the  electrical  properties  of  the  surface  and  the  local  geometry,  and 
another  factor  is  a  function  of  the  energy  of  the  small-scale  roughness  £(r');  refer  to 
Eq.  (20).  Thus,  if  the  illuminated  area  is  of  approximately  one  “patch”  (i.e.,  narrow  pulse 
and  narrow  antenna  beam  esse),  sea-clutter  statistics  may  be  modeled  by  a  product  of  two 
random  variables: 


Z  =  XY ,  (21) 

where  X  is  a  random  variable  dependent  on  the  slopes  of  the  large-scale  roughness,  £'  and  Y 
is  another  random  variable  dependent  on  the  statistics  of  the  small-scale  roughness  £.  The 
variables  X  and  Y  here  are  taken  to  be  mutually  independent. 

Accordingly,  the  probability  density  functions  (pdfs)  of  the  random  variables  X  and  Y 
should  be  of  the  forms 


p(*)=p«-') 


dT 

dX 


and 


p(Y)  =  e'Y, 


(22) 


(23) 


where  the  mean  value  of  Z  has  been  included  in  X. 

The  reason  for  selecting  these  two  pdfs  should  become  more  obvious  as  we  proceed. 
Equation  (22)  follows  directly  from  the  law  of  transformation  of  pdfs,  and  Eq.  (23)  is  the 
pdf  of  the  energy  spectral  density  of  the  small-scale  roughness,  which  is  assumed  to  be 
Gaussian  in  amplitude.  Therefore,  the  pdf  of  the  random  variable  Z,  according  to  elementary 
probability  theory,  should  be  given  by  (16) 

OO  J 

P(Z)  =  J  ~p(X)p  (r=|W.  (24) 

0 

For  the  case  involving  illumination  of  several  “patches”  at  one  time  (wide-pulse  and 
wide-antenna  beams),  sea-clutter  statistics  may  be  represented  by  the  sum  random  variable. 

Z(=Z1+Z2  +  ...  +  Z„.  (2B) 

Clearly,  the  conditional  density  p(V/H0)  used  in  calculating  the  probability  of  false 
alarm  in  Eq.  (2)  is  Eq.  (24)  for  the  narrow-pulse  and  narrow-antenna  beam  or  the  pdf  of  the 
random  variable  in  Eq.  (25)  for  the  wide  case. 

A  detailed  development  follows  for  the  simpler  case  of  single  “patch”  illumination. 

The  random  variable  X  in  Eq.  (21)  can  readily  be  identified  with  the  normalised  radar  cross 
section  obtained  for  slightly  rough  scattering  theory  (7),  thus  (to  first  order),  for  horizontal 
polarization 
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X(di)IJ  =  4 irft4  cot4  0,-|a2  cos2  6TXi  +  sin2  6[e(l  +  a2) - a2J  I’ll  M  I2 • 
W(2ka,2ky  sin  6), 

and  (to  first  order)  for  vertical  polarization 

X(&i)v  =  4irfe4  cot4  0,|a2  cos2  5[e(l  +  ctf)-ctf]Tu  ,,  +  sin2  67^  I2  • 


W(2ka,2ky  sin  5), 


where 


jj  =  cos  dj  “  cos  (8  +  }]/)  cos  5,  ctj  =  (1  -  yf)112,  a  =  sin  (0  +  il/),y  =  (1  -a2)1**, 

- C  —  - —  y<  =  ^ 

(Ti+T')2'"11  (E7,*-rV7  '  • 

6  is  the  complex  dielectric  constant  of  the  surface,  W{K x ,  Ky )  is  the  two-dimensional  spec¬ 
trum  of  the  small-scale  roughness  of  the  “patch”  (i.e.. 


o O 

^  =  (1/4)  II  W(KX  JCy  )  dKx  dKy  \ 


which  is  taken  to  be  constant  on  each  “patch,”  and  ^  and  6  are  the  “tilt”  angles  with  respect 
to  the  horizontal  plane,  parallel  and  perpendicular  to  the  plane  of  incidence,  respectively. 

To  have  a  unique  relationship  between  X  and  slope  f  \  side  tots  cannot  be  included  in 
this  formulation.  Thus,  in  this  statistical  model,  only  tilts  parallel  to  the  [Jane  of  incidence 
are  to  be  included.  Letting  5  =  0  in  Eqs.  (26)  and  (27),  gives 


X(9,)„  =  4irfe4  cos4  0;  -  - — ^7=== 
v  l  H  1  (c^  0,  +v/e-sm2  8,) 


W(2k  sin  8it  0) 


X(6;)v  -  4jrft4  cos4  Of  * *  W(2k  sin  8jt  0).  (29) 

(e  cos  8j  +  >/ e  -  sin4  0,-)4 

The  slopes  of  the  undulating  surface  are  Gaussian  distributed, 

(HO) 

and  since  tan  \p  =  the  pdf  of  X  can  be  found  from  Eq.  (22)  using  Eq.  (30)  and  Eqs.  (28) 
or  (29).  Let  W(KX,KY)  “  6  •  1(T3E“4  and  initially  take  the  case  of  a  perfectly  conducting 
sea,  (where,  e  m  Then,  Eqs.  (28)  and  (29)  become 
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X(dj)H  =  C  cot4  0, 


(31) 


and 


X(di)v  -  C  (cot2  9t  +  2)2, 


(32) 


where  C  *  (3ir/2)  •  lO-3. 

After  some  straightforward  algebra,  for  a  perfectly  conducting  sea  the  pdfs  for  the 
random  variables  X  are 


and 


_  1  /  A-tanr  \  2 
aec2  r  e  2S3\1+Atan T  j 

P(X/f)0C4 CS^EF  A3(l  +A  tan  r)2 


(33) 


P(*y)  a 


_  1  1 

f  B-tanr  \2 

e  M*1 

^1+Btanry 

4CS^2if  '  B(B2  +  2)(1  +  B  tan  t)2  ’ 


(34) 


where  r  *  90°  -  6  (grazing  angle),  A  =  and  B  *  [(Xy/C)1!2  -  2] 1^2.  The  con¬ 

stants  of  proportionality  in  Eqs.  (33)  and  (34)  are  obtained  from  the  normalization  of  the 
pdf. 


In  Figures  2  through  5,  Eq.  (33)  has  been  plotted  for  various  grazing  angles  and  rms 
slopes  of  the  undulating  surface.  As  observed  the  pdf  of  XH  is  nearly  Gaussian  for  small  S 
and  tends  to  the  exponential  distribution  for  large  S.  For  small  grazing  angles  the  tail  of  the 
pdf  increases  with  large  S,  a  characteristic  of  the  lognormal  distribution;  so  seemingly  we 
have  found  that  the  lognormal  properties  of  sea  clutter  may  be  explained  in  terms  of  the 
tilting  of  th"  slightly  rough  patches  by  the  large-scale  roughness. 


The  probability  distribution  corresponding  to  Eq.  (33)  is 

nr>JW«i[»r(sgJ-)  ■  m 

where  the  upper  sign  applies  for  tan  ^  >  0  ana  the  lower  sign  applies  for  tan  ^  <  0,  tan  ^  * 
A  -  tan  t/1  +  A  tan  r  and 


*/<-) 


dt 


Equation  (35)  has  been  plotted  in  Fig.  6  for  r  *  40°  and  for  several  values  of  S.  Figure 
7  gives  same  curves  plotted  on  normal  probability  paper,  where  lognormal  distributions  plot 
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Rif.  2— Density  function  of  Xjj  tar  perfectly  cooduetinf  eaa  and  40- 


as  straight  lines.  These  curves  indicate  that  the  distribution  of  XH  has  lognormal  character¬ 
istics,  which  is  most  evident  for  small  S,  and  for  large  S  the  tail  of  the  distribution  increases 
with  increasing  5. 

The  pdf  and  the  distribution  function  for  Xy  tor  a  perfectly  conducting  aea  are  not 
too  interesting  because  Eq.  (32)  approaches  a  constant  toward  the  grazing  incidence.  Thus, 
to  obtain  realistic  results  for  the  distributions  of  Xy  the  more  exact  expression  Eq.  (29), 
which  applies  for  a  sea  surface  of  finite  conductivity,  must  be  used.  The  pdf  and  the 
probability  distribution  for  X^  will  not  change  drastically  for  the  case  of  finite,  but  large 
conductivity. 

From  the  functional  dependence  of  Eqs.  (28)  and  (29)  on  the  angle  of  incidence  and 
tilt  angle,  the  inverse  of  those  expressions  can  be  approximated  by  the  aeries 

N 

tan(90°  -  0()  -  £  aiX*'N , 

£-1 


(36) 
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XH 


Fig.  3— Density  (unction  of  Xfj  for  perfectly  conducting  sea  and  30- 
degree  grazing  angle 


(38) 
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Pi*.  4— Density  function  of  Xfj  for  pwrfoctly  conducting  tea  and  20- 
d«grM  grazing  angle 


where  again  the  upper  sign  applies  for  tan  i //  >  0  and  the  lower  sign  applies  for  tan  \l/  <0. 
In  this  case 


tan  ^ 


N 


^  a%X*/N  -  tan 
£-1 


T 


N 

1  +  tan  r  ^  a^X^/N , 
£-1 


and  the  factor  of  proportionality  in  Eqs.  (37)  and  (38),  as  before,  is  obtained  from  the 
normalization  condition  P(X  >  0)  •  1. 

A  more  careful  investigation  of  Eqs.  (37)  and  (38)  show  that  the  distribution  of  X 
should  tend  toward  the  Gaussian  distribution  when  the  last  few  coefficients  dominate  (this 
is  the  case  for  vertical  polarization).  On  the  other  hand  if  the  first  few  coefficients  dominate 
in  Eq.  (36),  the  distribution  of  X  should  tend  toward  the  lognormal  distribution  (this  is  the 
case  for  horizontal  polarization.  Thus,  this  model  predicts  that  sea  clutter  for  horizontal 
polarization  should  tend  toward  the  lognormal  distribution. 


i 
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Fig.  5— Density  function  of  Xf]  for  perfectly  conducting  eee  and  10- 

degree  gracing  angle 


Of  course,  the  statistics  of  sea  clutter  in  final  analysis  is  related  to  the  distribution  of 
the  random  variable  Z  which  also  depends  on  the  distribution  of  the  other  random  variable 
Y,  which  in  our  model  is  always  an  exponential  distribution.  Thus,  any  lognormal  char¬ 
acteristics  of  sea-clutter  distribution  must  come  from  the  distribution  of  the  random  variable 
X.  In  general  closed-form  expressions  cannot  be  obtained  for  the  pdf  and  distribution  func¬ 
tion  of  Z. 

However,  some  limiting  properties  can  be  derived  for  the  distribution  of  Z.  For  example, 
take  the  cue  of  a  glassy  aea  (no  undulating  function  is  present,  S  -  0).  Clearly  the  pdf  of 
X  will  be  a  delta  function: 

(39) 


and 


(40) 


where  7  is  the  mean  value  of  X  and  also  of  Z,  and  Z  is  exponentially  distributed  for  this  cat*. 


Of  EXCEEONG  THE  ABSCISSA  *  KX>  PROBABILITY  OF  l-XCEEOMG  T* 
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For  increasing  value  of  5  (rougher  seas),  p(X)  should  develop  a  longer  tail  (see  Figs. 
2-5);  thus,  p(Z)  will  also  have  a  longer  tail,  which  should  make  the  distribution  of  Z  tend 
toward  the  lognormal  distribution  (in  particular  for  horizontal  polarization). 

The  moments  of  Eqs.  (33),  (34),  and  (37)  do  not  exist  because  we  have  used  a  K~ 4 
spectrum  which  has  a  singularity  at  K  ■  0.  This  in  principle  can  be  corrected  by  cutting  off 
the  spectrum  before  K  m  0.  A  shadowing  function  may  also  be  included  in  Eq.  (21 )  but 
has  not  been  attempted  here. 


IDENTIFICATION  OF  THE  STATISTICS  OF  SEA  CLUTTER  OBTAINED 
WITH  THE  NRI-4FR  SYSTEM 

Quantitative  analysis  of  sea  clutter  statistics  is  necessary  for  further  development  and 
improvement  of  the  statistical  model  of  sea  clutter.  Thus,  analysis  has  been  performed  on 
•es -clutter  data  taken  with  the  4FR  system  which  generally  uses  pulse  widths  between  0.1 
and  0.5  psecond  and  whose  antenna  beam  widths  are  greater  than  5  degrees.  Accordingly 
the  illuminated  area  in  this  case  contains  many  “patches”  and  no  direct  comparisons  will  be 
possible  with  the  predictions  of  the  statistical  model,  but  any  statistical  information  obtained 
for  sea  clutter  will  be  helpful  in  assessing  or  updating  the  statistical  model. 

A  detailed  description  of  the  4FR  system  has  already  been  given  by  Guinard  (8);  thus, 
here  we  should  only  indicate  that  the  samples  of  sea  clutter  used  in  this  investigation  are 
samples  of  backscattered  power  (or  radar  cross  section  of  the  sea)  which  have  been  collected 
by  a  logarithmic  receiver  with  a  dynamic  range  of  over  45  decibels.  The  output  of  the  re¬ 
ceiver  is  digitalized  by  means  of  a  30-nanosecond  gate  to  seven-bit  accuracy.  The  sea  condi¬ 
tions  and  radar  parameters  of  the  measurements  are  summarized  in  Table  1.  The  experiment 
in  the  North  Atlantic  during  February  1969  was  performed  in  the  neighborhood  of  Ocean 
Stations  India  and  Juliet. 

The  Distribution  Function 

The  samples  of  sea  clutter,  for  depression  angles  of  5  to  30  degrees,  taken  with  the 
4FR  system  have  been  distributed  and  checked  with  the  Kolmogorov-Smimov  test  (17)  if 
they  come  from  an  exponential  population  (Rayleigh  envelope)  or  a  lognormal  population. 

In  fitting  the  exponential  distribution,  the  mean  value  of  the  exponential  distribution  was 
adjusted  for  a  minimum-maximum  deviation  between  the  sample  and  the  population  dis¬ 
tribution.  In  fitting  the  lognormal  distribution,  the  variance  and  the  median  values  were 
taken  to  be  those  of  the  sample. 

‘Die  number  cf  independent  samples  of  sea  ciutter  was  estimated  by  means  of  a  stand¬ 
ard  run-test  (18)  in  an  application  of  the  Wald-Wolfowitz  test  (19).  The  procedure  used 
here  is  identical  to  that  described  by  Schmidt  (20)  except  that  in  our  investigation  sets  of 
1024  samples  are  used  instead  of  2048.  The  run-test,  using  the  95%  level  of  significance, 
shows  that  1  out  of  4  samples  is  independent  for  X  band,  1  out  of  8  or  9  for  C  band,  1  out 
of  35  for  L  band,  and  1  out  of  120  to  130  is  independent  for  P  band.  Of  course,  these  are 
average  numbers;  in  the  analysis  the  exact  number  of  independent  samples  was  used  for  each 
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Table  I 

Radar  and  Sea  Conditions 


Date 

Wind  Speed 
(m/sec) 

Wave  Height 
(m) 

Location 

Altitude 

(m) 

Indicated 
Air  Speed 
(m/sec) 

Pulse 

Rep. 

Rate 

(PPa) 

Pulse  Width 
(nsec) 

2/11/69 

24.5 

8.0 

OWS  “I” 

480 

1 

683 

0.6 

2/13/69 

18 

OWS  “J” 

460 

■ 

603 

0.5 

2/14/69 

20 

7.5 

OWS  “J” 

460 

603 

0.6 

2/17/69 

2.5 

4.6 

OWS  “I” 

410 

103 

603 

0.5 

2/18/69 

11 

3.0 

OWS  “J” 

460 

103 

603 

0.6 

1/23/70 

6.2 

3.7 

Bermuda 

180 

103-115 

683 

0.6 

1/26/70 

7.5 

1.5 

Bermuda 

750 

89-108 

683 

0.5 

1/27/70 

8.2 

1.8 

Bermuda 

180 

87-108 

683 

0.5 

1/27/70 

8.2 

1.8 

Bermuda 

_ 

610 

87-110 

683 

0.5 

cumulative  distribution  and  in  general  the  number  of  independent  samples  is  proportional  to 
the  radar  frequency.  No  specific  trends  were  found  with  polarization;  the  motion  of  the 
aircraft  is  expected  to  decorrelate  the  samples  of  sea  clutter  faster.  Thus,  for  sea  clutter 
taken  with  a  stationary  radar  the  number  of  independent  samples  should  be  smaller. 

The  Kolmogorov-Smimov  test,  either  at  the  99%  or  at  the  80%  level  of  significance, 
shows  that  for  calm  seas  and  small  sample  sizes  the  cumulative  distributions  of  sea  clutter 
are  acceptable  as  exponential  or  lognormal.  However,  for  large  sample  sizes  (data  up  to  30 
seconds  in  this  case)  the  maximum  deviation  between  the  distributions  becomes  larger  than 
the  acceptable  limit  and  should  be  rejected  as  belonging  to  the  exponential  or  the  log¬ 
normal  family,  with  the  smaller  maximum  deviations  obtained  against  the  exponential  dis¬ 
tribution  in  most  cases. 

The  maximum  deviations  for  P-  and  L-band  data  were  in  general  smaller  than  the  maxi¬ 
mum  deviations  for  C-  and  Jf-band  data  when  compared  with  the  exponential  and  the  log¬ 
normal  distributions.  A  trend  occurs  with  polarization  which  may  be  more  obvious  in  terms 
of  the  central  moments,  which  are  investigated  in  the  next  section. 

In  Figs.  8a-8x,  some  typical  maximum  deviations  illustrate  the  above  conclusions.  In 
these  figures,  only  the  comparison  with  the  exponential  distribution  is  shown.  The  maximum 
deviations  in  comparison  with  the  lognormal  distribution  are  not  too  different  than  those 
shown,  except  that,  in  general,  they  tend  to  be  a  little  greater. 

The  outcome  of  the  Kolmogorov -Smirnov  test,  although  conclusive,  is  not  very  satis¬ 
fying  since  the  statistics  of  sea  clutter  are  still  not  known. 
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The  Central  Momenta 

rrhe  Kolmogorov-Smimov  test  has  demonstrated  that  sea  clutter  taken  with  the  4Fft 
system  is  not  xponential  nor  lognormal  distributed.  To  obtain  a  quantitative  estimate  of 
these  statistics,  the  fust  five  central  moments  of  sea  clutter  should  be  computed  from  the 
data,  in  decibels,  to  -.uud  additional  errors  in  scaling  of  the  data. 

The  central  moments  Pj ,  M 2.  — .  M„  of  a  random  variable  which  is  exponentially  or  log¬ 
normal  distributed  are  well  known  and  can  be  used  for  comparison.  The  central  moments 
for  the  logarithm  of  a  random  variable  which  is  exponentially  distributed  are  related  to  the 
poligamma  function  (21) 


(where  F(x  +  1)  is  the  gamma  function).  The  central  moments  are: 

p2  =  ^'  (0)  >»  1.64493  > 

=  (0)=-  -  2.40411  \ 

p4  =  ^'•1  (0)  =  6.49393 

p5  =  i//IV  (0j  = -24.88627  , 


(41) 


(42) 


and  \li(0 )  =  -0.5772157  is  related  to  the  difference  hetw*en  the  natural  logarithm  of  the 
mean  value  of  the  exponential  distribution  and  the  mean  value  of  the  natural  logarithm  of 
the  random  variable.  The  numerical  values  of  the  poligamma  function  have  been  obtained 
from  Davis  (22). 


The  central  moments  of  the  logarithm  of  a  random  variable  which  has  a  lognormal  dis¬ 
tribution  are  those  of  the  Gaussian  distribution  (14) 

p„  =  l-3-5...(n-l)pa»/2,  (43) 

where  n  >  2  and  is  even.  If  n  is  odd,  p„  *  0. 

In  Figs.  9  and  10  the  puf  and  the  distribution  function  of  an  exponential  and  a  log¬ 
normal  random  variable  have  been  sketched  for  comparison. 

Figures  11  through  20  illustrate  the  (median,  mean)  difference  and  the  first  five 
central  moments  of  sea  clutter  (in  decibels)  taken  with  the  NRL-4FR  system  for  moderate 
and  rough  sea  conditions.  These  figures  also  show  the  corresponding  values  for  the  exponen¬ 
tial  and  lognormal  distribution.  As  observed  the  central  momenta  of  sea  clutter  are  in  most 
cases  intermediate  between  those  of  the  exponential  (Rayleigh  envelope)  and  those  of  the 
lognormal  distribution,  except  for  the  fourth  central  moment. 
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(f)  Horizontal  polarization  of  L  band  on  Feb. 
11,1969 


(e)  Horizontal  polarization  of  P  band  on  Feb. 
11, 1969 


Fig.  8— Maximum  deviation  of  cumulative  diatri button  of  tea  clutter,  taken  with  the  4PR  tj (tern,  from  the 
exponential  diitribution  (Rayleigh  envelope),  with  the  96%  and  60%  leveia  of  rignifVance  from  the 
Kolmogorov- Smirnov  teat 


(h)  Horizontal  polarization  of  X  band  on  Feb. 

11. 1969 


K 


UM  DEVIATION  OF  CUMULATIVE  DCSTTWU 
FROM  THE  EXPONENTIAL  DRTRfiUTlON 


2 


99% 
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(p)  Horizontal  polarization  of  X  band  on  Jan. 
26, 1970 


(r)  Vertical  polarization  of  L  band  on  Jan.  26, 
1970 


(q)  Vertical  polarization  of  P  band  on  Jan.  26, 
1970 


Fig.  8— Maximum  deviation  of  cumulative  distribution  of  tea  clutter,  taken  with  die  4PR  system,  from  the 
exponential  distribution  (Rayleigh  envelope),  with  the  99%  and  80%  levels  of  significance  from  the 
Kolmogorov-Smirnov  test 


IV 


UM  OCVIATON  OF  CUMULATIVE  OtSTfW. 
FROM  T^C  EXPONENTIAL  QSTRBUTION 


NRL 


(v)  Horizontal  polarization  of  L  band  on  Jan. 
28,  1970 


(z)  Horizontal  polarization  of  X  band  on  Jan. 
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Fig.  9— Compariton  of  the  exponential  (Rayleigh  envelope)  and 
the  lognormal  denaity  functions 


Fig.  10-  Compariion  of  tliu  exponential  (Rayleigh  envelope)  and 
the  lognormal  distribution  function* 


The  central  momenta  for  vertical  polarization  are  nearer  to  those  for  the  exponential 
distribution,  which  agrees  with  the  prediction  of  the  model  for  the  statistics  of  sea  clutter. 

It  would  seem  that  the  wide  spread  in  the  magnitude  of  the  central  moments  may  be  an  indi¬ 
cation  that  in  general  sea  clutter  is  not  stationary.  Temporal  variation  of  sea  clutter  has  been 
investigated  by  Schmidt  (20). 

With  the  central  moments  obtained  for  sea  clutter,  the  true  pdf  and  distribution  func¬ 
tions  should  be  reconstructable  with  no  difficulty. 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

A  model  /or  the  statistics  of  sea  clutter  has  been  developed  using  acattering  theory 
and  the  composite  surface-scattering  model,  which  uses  the  assumption  that  the  rough 
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R|.  11— Median,  mean  difference  of  aea  clatter  taken  with  the 
4FR  system  for  moderate  aea*.  The  wind  speed  waa  5  to  10 
meters /second.  Solid  point*  repreeent  vertical  polarization  and 
open  points,  horizon taJ  polarization. 
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Fig.  12— Median,  mean  difference  of  aea  clatter  taken  with  the 
4FR  system  for  rough  seas.  The  wind  speed  was  15  to  20 
meters/second.  Solid  points  represent  vertical  polarization  and 
open  points,  horizontal  polarisation 
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Fig.  13— Variance  of  tea  clutter  taken  with  the  ■»FK  ryaiem  for 
moderate  aeaa.  The  wind  apeed  «u  5  to  10  metera/aecond.  The 
aolid  pointa  repreaent  vertical  polarization  and  the  opei>  point*, 
horizontal  polarization. 


KAMft  raeaucMCir  (ucahckt;) 


Fig.  14— Variance  of  aaa  chatter  taken  with  the  4FR  ajratean  for 
rough  aaaa.  TV*  wind  apeed  waa  16  to  SO  owtan/aaoosd.  TV* 
aolid  point*  rapr.aent  aartkai  polarization  and  open  point*, 
hortaootol  polamahoa. 
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RADAR  FREQUENCY  (GISAHERTZ) 

Fig.  17— Fourth  central  moment  of  sea  clutter  taken  with  the 
4FR  system  for  moderate  seas.  The  wind  speed  was  5  to  10 
meters/second.  The  solid  points  represent  vertical  polarization 
and  the  open  points,  horizontal  polarization. 


RADAR  FREQUENCY  i GIGAHERTZ) 

Fig.  18— Fourth  central  moment  of  sea  clutter  taken  with  the 
4FR  system  for  rough  seas.  The  wind  speed  was  15  to  20 
metei/second.  The  sc  lid  points  represent  vertical  polarization 
and  the  open  points,  horizontal  polarization. 
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surface  “locally”  is  slightly  rough.  For  the  sea  the  “local”-slightly  rough  assumption  is 
quite  a  reasonable  approximation  at  most  conventional  radar  frequencies. 

The  case  of  single  “patch”  illumination  is  treated  in  detail,  with  the  model  predicting 
the  statistics  of  sea  ciutter  to  be  exponential  (Rayleigh  envelope)  for  a  glassy  sea  (no  large- 
scale  roughness  present),  and  for  increasing  roughness  the  distribution  should  develop  a 
longer  tail  similar  to  that  for  a  lognormal  distribution,  in  particular  for  horizontal  polariza¬ 
tion.  Obviously,  a  mechanism  has  been  found  that  yields  sea  clutter  with  a  distribution  of 
lognormal  character. 

Beckmann  (23)  has  shown  that  by  using  physical  optics  for  a  perfectly  conducting 
rough  surface  the  scattered  fields  should  be  Rayleigh  distributed  everywhere  except  near  the 
specular  direction,  and  for  very  rough  surfaces  the  scattered  fields  should  be  Rayleigh  dis¬ 
tributed  even  in  the  specular  direction.  Thus,  the  non-Rayleigh  statistics  are  obviously  due 
to  the  scattering  properties  of  the  slightly  rough  “patches’*  and  the  tilting  of  the  “patches” 
by  the  large-scale  roughness. 

As  future  input  for  updating  the  model,  the  statistics  of  sea  clutter  taken  with  the  4FR 
system  have  been  analyzed  by  the  Kolmogorov-Smimov  test  and  by  a  computation  of  the 
central  moments  of  the  distribution.  The  results  indicate  that  in  general  sea  clutter  is  not 
exponentially  (Rayleigh  envelope)  nor  lognormally  distributed,  and  these  distributions  may 
only  be  the  limiting  distributions  of  sea  clutter.  The  large  spread  of  the  central  moments 
may  be  an  indication  that  sea  clutter  is  not  stationary;  therefore,  the  optimum  radar  detec¬ 
tor  must  be  of  an  adaptive  nature. 

Some  other  important  results  obtained  in  the  empirical  identification  of  sea  clutter  are: 

1.  For  large  sample  sizes  (greater  than  about  1000  independent  samples)  sea  clutter 
is  not  exponential  nor  lognormal, 

2.  For  small  sample  sizes  (less  than  200  independent  samples)  sea-clutter  statistics 
may  be  approximated  by  both  the  exponential  and  the  lognormal  distributions, 

3.  Sea  clutter  for  vertical  polarization,  in  general,  is  more  exponential  than  sea  clutter 
for  horizontal  polarization, 

4.  Sea  clutter  for  P  and  L  band,  in  general,  is  more  exponential  than  sea  clutter  for  C 
and  X  band, 

5.  The  number  of  independent  samples  is  roughly  proportional  to  radar  frequency, 

6.  Sea  clutter  for  calm  seas  is  more  exponential  than  sea  clutter  for  rougher  seas. 

In  conclusion,  it  is  possible  to  say  that  the  model  developed  for  the  statistics  of  sea 
clutter  and  the  empirical  identification  of  data  taken  with  the  4FR  system  indicate  that  the 
distribution  of  sea  ciutter  is  intermediate  between  the  exponential  (Rayleigh  envelope)  and 
the  lognormal  distribution.  Accordingly,  the  expected  probability  of  false  alarm  for  radars 
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Fig.  19— Fifth  central  moment  of  sea  clutter  taken  with  the  4FR  ! 

system  for  moderate  teas.  The  wind  speed  wac  5  to  10  meters/ 
second.  The  solid  points  represent  vertical  polarisation  and  the 
open  points,  horizontal  polarization. 
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Fig.  20— Fifth  central  moment  of  sea  clutter,  taken  with  the  4FR 
system,  for  rough  seas.  The  wind  speed  was  15  to  20  meters/ 

second.  The  solid  points  represent  vertical  polarization  and  the  ! 

open  points,  horizontal  polarization. 
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operating  in  the  sea  environment  should  be  larger  than  that  predicted  for  an  exponentially 
distributed  clutter,  in  particular  for  horizontal  polarization  and  toward  grazing  incidence. 
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